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Abstract 

Interchain hopping of solitons and polarons in polyacetylene is studied by nu- 
merical simulation of their motion under an electric field. Use is made of the 
Su-Schrieffer-Heeger model supplemented with intrachain electron-electron in- 
teractions and dopant potentials. We find that charged solitons can hop to 
the opposite chain by forming bound pairs (bipolarons). For the case of po- 
larons also, hopping in a pair is more favorable than single polaron hopping. 
Interchain hopping of a polaron pair is more efficient than that of a soliton 
pair. 

I. INTRODUCTION 

The mechanism of charge hopping among chains is one of the intricate unsettled ques- 
tions related to understanding of the electronic properties of conducting polymers. |3J It is 
expected that charged solitons or polarons play a major role in conduction processes in poly- 
acetylene. Motion of such nonlinear excitations can contribute directly to the conduction. 
PU Polyacetylene, however, has a complex structure composed of polymer chains with finite 
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lengths and dopant ions (Fig. 1). A difficulty in the charged soliton model is that a single 
soliton cannot hop to neighboring chains, since there exists a large activation barrier for 
inversion of the phase of bond alternation over a large region. The difficulty can be removed 
by consideration of pairs of solitons, i.e. bipolarons. Interchain hopping of polarons, which is 
another conduction mechanism, is considered to be much easier. It is crucial to consider the 
dynamics of the nonlinear excitations by taking into account the three-dimensional effects, 
such as dopant potentials and interchain interactions. 

Silva and Terai Q investigated the effect of the interactions between two chains of poly- 
acetylene on the dynamics of solitons, and found that the motion of a soliton is suppressed 
by a confinement effect and that the amount of the charge transfer between the two chains 
via charged and neutral solitons is quite small. The effect of a dopant has been usually 
incorporated as the Coulomb potential for electrons. Yamashiro and coworkers pointed out 
that interchain transfers are enhanced by the presence of the valence orbitals of the dopant 
ions. U 

We study here interchain hopping of solitons and polarons in polyacetylene by numerical 
simulation of their motion under an electric field. We consider two neighboring chains that 
are coupled with each other through dopant ions as shown in Fig. 1(b). One of these chains 
contains a pair of solitons or polarons in the initial state. The effects of a dopant are twofold: 
the Coulomb potential and enhancement of interchain transfer around the dopant. In our 
model, the interchain transfer is assumed to occur in a small region around the dopant, and 
is neglected in the other regions, for simplicity. 

The effect of chain ends on the motion of the nonlinear excitations is also important. 
For discussion of this effect the dopant is located near the end of one of the chains, chain 
a in Fig. 1(b). The other chain (chain b) is arranged so that the dopant is located near 
its midpoint. Thus we choose the structure shown in Fig. 1(b), and it seems reasonable to 
regard this structure as a typical example of part of the three-dimensional complex structure 
shown in Fig. 1(a). 
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II. MODEL AND NUMERICAL METHOD 



The model Hamiltonian is given by 

H = H a + H h + H e _ c + H mt + H imp . (2.1) 

The first two terms, describing the two chains a and b, are the Su-Schrieffer-Heeger model 
H modified to include the electric field. 

H\ = -XX*o - oty x ,i) [e 17A 4,i<7 C A,i+i CT + h.c] 

ia 

+ f YAi + yE<< - *W - «A,l)> (2-2) 
i j 

where A = a, b and — u x,i+i — u x,i- Here, to is the hopping integral between the nearest 
neighbor sites, a the electron phonon coupling constant, ux,i the displacement of the z'th 
CH unit on the chain A, c\^ a the annihilation operator of an electron at the z'th site with 
spin cr on the chain A, K the spring constant, and M the CH mass. The electric field is 
applied in a direction parallel to the chain axis, being given by the time dependent vector 
potential, E = —A/c. The parameter 7 is defined by 7 = ea/hc, where e is the absolute 
value of the electron charge, a the lattice constant, and c the light velocity. We adopt the 
open boundary condition for each chain. The last term in eq. Q2.2|) is added to keep the 



chain length constant. The third term in eq. (|2.1|) describes the short-range intrachain 
electron-electron interactions. 

H c ^ e = UJ2 n x,i\ n \H + V J2 n x,in\, i+ i, (2.3) 

X,i X,i 

with nx t ia = c Xia C\ t i(j and ri\^ = J2a n x,ia- The parameters U and V are the on-site and 
nearest neighbor Coulomb repulsion, respectively. The two chains are of the same size. 
Chain b is shifted by z'q sites with respect to chain a. The interaction between two chains is 
represented by 
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where p and q denote the first and last site in the interaction region around the dopant. The 
last term of eq.(|2.1|) represents the dopant potential, which is approximated by a square-well 
potential, 

<? 

H mi p = V[ mp ( C a,icr C a,i(T + c b,i-i a C ^,i-io^) ' (2-5) 

<J,i=p 

We use the following values of parameters: |||| — 2.5eV, t± = 0.1t , K = 21eV/A 2 , 
a = 4.1eV/A, a = 1.22A, V imp = -0.1t , U = t , V = t /2 and T = 5.13eV. The chain 
sizes are iV a = iVb = 120 and the total number of electrons is N e = 242, corresponding to 
two excess electrons in the system. The center of the dopant is located at the 105th site 
in the coordinate of chain a. The endpoints of the interaction region are set to be p = 100 
and q = 110. The shift of chain b with respect to chain a is %q = 46, so that the dopant is 
located near the center of chain b. 

The initial states are prepared self-consistently in the Hartree-Fock approximation, in 
which chain a contains two solitons or polarons, one of which is pinned by the dopant 
and the other of which is free to move. The time evolution of the electronic wavefunction 
is determined by use of the time-dependent HF equation, and the motion of the lattice 
by use of an equation of motion. In solving the time dependent HF equation we use 
the method of fractal decomposition for exponential operators. [ 10|, 1 1 1 For the equation of 



motion the time differential equation is integrated with discretization of the time with an 
interval At. The interval At is chosen to be At = 0.0025cg>q 1 throughout this work. Here in 
our simulations the unit of time is the inverse of the bare optical phonon frequency (uq = 



AK/M m 2.5 x lO^V 1 ) and that of the electric field is E = hu/ea (w 1.3 x 10 7 V/cm). 



III. RESULTS AND DISCUSSION 



A. Case of a pair of solitons 



In the following we show the results in terms of the smoothened quantities of the bond 
variable and the excess charge distribution, defined as 
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y x ,i = (-1)%a,-i - 2yx,i + 2/A,i+i)/4, (3.1) 
P\,i = (Px,i-i + 2p A ,i + PA,i+i)/4, (3.2) 

where the excess charge density is calculated from the time-dependent wavefunctions 

{*A*(M)}, 

/M*) = £V;u,(M)| 2 -l- (3.3) 

v 

The prime attached to the summation symbol denotes a sum over occupied states. 

Figure 2 shows a typical example of hopping of a pair of solitons from chain a to b in 
three-dimensional representation of y\i an d p\i as functions of time and space. The free 
soliton on chain a is accelerated by the field and collides with the pinned soliton, resulting 
in formation of a bipolaronic state. In chain b small oscillations are generated after the field 
is switched on. Their amplitudes increase upon the collision and sometimes the oscillations 
induce bipolaronic lattice distortion. In such a case the soliton pair can easily hop to the 
neighboring chain. After the hopping, one of the solitons is pinned by the dopant and the 
other begins to move toward the chain end. 

We have performed numerical simulations for various strengths of the field for changing 
of the incident velocity of the free soliton at the collision with the pinned soliton. In Fig. 3 
we display the time dependence of the excess charge density on chain a, which is defined by 

Pa(t)=EPa,(t). (3.4) 
i=l 

In the initial state p&(0) has almost 2|e|, since there exist a pair of charged solitons on chain 
a. A slight deviation from 2|e| is due to partial flow of charge from the pinned soliton to 
chain b. The deviation increases and fluctuates after the two solitons form a bipolaron. In 
the case of a weak field (i.e. low incident velocity at the collision), the soliton pair does 
not show any sign of hopping at least within the longest time of our simulations. For a 
stronger field (high incident velocity), p a (t) falls to about zero suddenly, indicating that the 
hopping of two electrons occurs simultaneously. There is a threshold velocity for the hop. 
This means that there is a potential barrier for the hopping. The time interval between 
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the collision and completion of the hop is not a simple function of the field strength (the 
incident velocity), but a rather chaotic one. The occurrence of the hop depends on the 
mutual motion of the soliton pair on chain a and the lattice oscillation on chain b. Namely, 
the hopping occurs under the conditions that the amplitude of the bipolaron becomes small 
on chain a and instantaneously becomes large on chain b. Therefore a slight difference in 
parameters results in a drastic change in the hopping time. Bipolaronic lattice distortion 
on the opposite chain is necessary for the hop to occur. The lattice fluctuations at finite 
temperature may affect the hopping, perhaps they make the hop easier. 

The location of the dopant affects the hopping probability. If the dopant is located near 
the midpoint of chain a, the probability of hopping of the soliton pair is very low. When 
the field is sufficiently strong, the pinned soliton is pushed out by its collision with the free 
soliton and begins to move within chain a. In the case that the dopant is located near the 
chain end, even if the pinned soliton is pushed out, it is immediately reflected at the chain 
end, resulting in the formation of a bipolaron again. Bipolaron formation is necessary for the 
hopping to occur. Therefore, the hopping probability is larger when the dopant is located 
near the chain end. 

B. Case of a pair of polarons 

For the case of a pair of polarons three-dimensional representations of y\i and p\i are 
shown in Fig. 4. The hopping process seems qualitatively different from that for the soliton 
pair. The polaron pinned by the dopant in the initial state hops back and forth between the 
two chains in the electric field. The bipolaronic state observed in the case of the soliton pair 
does not appear as an intermediate state during the hop. The intermediate state in this case 
is such that each chain contains one polaron around the dopant. After the hopping, one of 
the two polarons is pinned by the dopant over the two chains and the other one is moved on 
chain b by the field. For the case of a single polaron, hardly any charge transport over the 
two chains due to interchain hopping through the dopant takes place, because of trapping 



6 



by the dopant potential. In Fig. 5 we show the time dependence of the excess charge density 
on chain a p & {t) for several field strengths. Before the collision, p a (t) oscillates at around 
1.5|e| due to the back- and- forth motion of the polaron pinned over the two chains. For 
the case E = 0.01-Eo ; two polarons are pinned by the dopant after the collision and there 
exists one polaron per chain. Therefore amount of the excess charge density on chain a is 
|e| on average. This case does not contribute to the charge transport from chain a to chain 
b. In the stronger field cases the charge transport occurs by the polaron hopping. Then 
p a decreases by |e| on average. The potential barrier for the hopping of a polaron pair is 
smaller than that for the hopping of a soliton pair. In contrast to the case for a soliton pair, 
the polaron pair hops almost irrespective of the dopant location relative to the chain end. 
When the two polarons collide, the hopping occurs almost instantaneously. Therefore, the 
time of the hop varies monotonically with the field strength, not in a chaotic manner. 

It seems that the back-and-forth motion of the pinned polaron before the collision is 
induced mainly by the external field. For the parameter set used in our simulations, the 
energy difference between the state of the polaron being almost localized on one chain and 
the state of delocalization over the two chains is quite small. Therefore the oscillation can 
be induced by weak disturbance. This is considered to be one of the reasons for the small 
potential barrier for the polaron pair hopping. 

It is worth noting that the final states result in formation of a pair of solitons or a pair 
of polarons depending on whether the initial state is a pair of solitons or pair of polarons. 
For a single chain with no electron-electron interaction, it has been shown by use of the 



TLM model El that a two-polaron state is unstable against a two soliton-state. [13] This 



is supported by the results of a numerical study on the dynamics of a polaron pair: it was 



observed that two polarons dissociate into a pair of solitons after their collision. In our 
two-chain system, a two-polaron state, not a two-soliton state, appears after the collision 
(hopping) of a pair of polarons. One of the reasons for this may be that the delocalization 
of polarons over the two chains stabilizes the polaron pair state. 

In summary, we have numerically investigated interchain hopping of a pair of solitons 



and a pair of polarons. We found that charged solitons can hop to the opposite chain by 
forming a composite (bipolaron) in the electric field parallel to the chains. The existence of 
both the dopant and the chain end assists the hopping. The bipolaronic lattice distortion on 
the opposite chain is necessary for the hop to occur. There exists a finite potential barrier 
for the hopping. For the case of polarons also, hopping in a pair is more favorable than 
single polaron hopping. Interchain hopping of a polaron pair is more efficient than that of 
a soliton pair. 
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Fig. 1. Schematic illustration of (a) polymer chains with dopants and (b) the system studied 
here. 
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Fig. 2. Stereographic representation of the site (i) and time (t) dependences of (a) the bond 
variable and (b) the excess charge density for the case of a pair of solitons. The interchain 
interaction region is depicted as the region between two lines in the bottom i — t plane. 

Fig. 3. Time dependence of the excess charge density on chain a for a pair of solitons. 

Fig. 4. Stereographic representation of the site (i) and time (t) dependences of (a) the bond 
variable and (b) the excess charge density for the case of a pair of polarons. 

Fig. 5. Time dependence of the excess charge density on chain a for a pair of polarons. 
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